A phenomenological theory of zero-energy Andreev resonant states 
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A conceptual consideration is given to a zero-energy state (ZES) at the surface of unconventional 
superconductors. The reflection coefficients in normal-metal / superconductor (NS) junctions are 
calculated based on a phenomenological description of the reflection processes of a quasiparticle. 
The phenomenological theory reveals the importance of the sign change in the pair potential for 
the formation of the ZES. The ZES is observed as the zero-bias conductance peak (ZBCP) in 
the differential conductance of NS junctions. The split of the ZBCP due to broken time-reversal 
symmetry states is naturally understood in the present theory. We also discuss effects of external 
magnetic fields on the ZBCP. 

PACS numbers: 74.50.-|-r, 74.25.Fy,74.70.Tx 



I. INTRODUCTION 

Transport phenomena in unconventional supercon- 
ductors have attracted considerable interst in recent 
years because high- Tc superconductors may have the 
d wave pairing symmetryiSi^ The unconventional 
pairing symmetry causes the anisotropy in transport 
properties such as the electric conductance and the 
thermal conductivity.^*Si In normal-metal / high- Tc 
superconductor junctions, for instance, the shape of 
the differential conductance reflects the density of 
states when the a axis of high- Tc materials is per- 
pendicular to the junction interface. When a axis 
deviates from the interface normal, on the other hand, 
the conductance shows a large peak at the zero bias- 
voltagei'^/.8,9,iQ,ii,i2,i3,i4,i5,i6,i7,i8 guch anisotropy in 

the conductance is now explained by the formation of a 
zero-energy state (ZES)^^ at the interface of junctions. 
Since the ZES appears just on the Fermi energy, it dras- 
tically affects transport properties through the interface 
of unconventional superconductor junctions. The low- 
temperature anomaly of the Josephson current between 
the two unconventional superconductors is explained 
in terms of the resonant tunneling of Cooper pairs via 
the ZES;20-li^_2.23,24 ^^^^ ^ considerable number of 

studies have been made on the ZES itself and related 
phenomena of transport properties in both spin-singlet 
and spin-triplet unconventional superconductor junc- 

7.8.25.26.27.28.29.30.31.32.33.34.35.36.37.38.39.40.41.42.43.44.45. 

The conductance in normal-metal / superconductor 
(NS) junctions is calculated from the normal and the An- 
dreev reflection^ coefficients which are obtained by solv- 
ing the Bogoliubov-de Gennes (BdG) equation-^ under 
appropriate boundary conditions at the junction inter- 
face. Consequently we easily find the zero-bias conduc- 
tance peak (ZBCP) in NS junctions of high- Tc supercon- 



ductorsi^ Although the algebra itself is straightforward, 
it is not easy to understand the physics behind the calcu- 
lation. In a previous paper, we briefly discussed reasons 
for the appearance of the ZBCP by a phenomenological 
argument. The phenomenological analysis has several 
advantages. For instance, it shows the importance of 
the unconventional pairing symmetry for the formation 
of the ZES without directly solving the BdG equation. 
Moreover we easily understand that the ZES is a result 
of the interference effect of a quasiparticle. The applica- 
bility of the analysis in the previous paper, however, is 
very limited because of its simplicity. 

In this paper, we reconstruct the phenomenological 
theory of the Andreev reflection to meet the mathemati- 
cal accuracy. We calculate the reflection coefficients of an 
electronlike quasiparticle incident from a normal metal 
into a NS interface. Near the junction interface, a quasi- 
particle suffers two kinds of reflection: (i) the normal re- 
flection by the barrier potential at the NS interface and 
(ii) the Andreev reflection by the pair potential in the 
superconductor. In the present theory, we consider the 
two reflections separately to calculate the transport coef- 
flcients. As a consequence, the Andreev reflection coeffi- 
cient is decomposed into a series expansion with respect 
to the normal reflection probability of NS junctions. The 
expression of the Andreev reflection probability enables 
us to understand the importance of the unconventional 
^airin^ symmetry for the formation of the ZES. In uncon- 
' uLiiliuliul superconductors, the pair potential in the elec- 
tron branch (A4.) differs from that in the hole branch 
(A_). The Andreev reflection probability at the zero- 
energy is expressed as the summation of the alternating 
series when A_|_ and A_ have the same sign with each 
other. In this case, the zero-bias conductance becomes a 
small value proportional to I^Afl*, where |tArp is the nor- 
mal transmission probability of junctions. On the other 
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hand when A+A_ < 0, all the expansion series have the 
same sign and the conductance has a large peak at the 
zero-bias. The phenomenological theory can be applied 
to superconductors with a broken time-reversal symme- 
try state (BTRSS)54,55,56,57,58,59,60,61,62,63,64,65,66,67 ^nd 

NS junctions under external magnetic fields i^^*S2iS2iZi 

This paper is organized as follows. In Sec. II, the An- 
dreev and the normal reflection coefficients are derived 
from a phenomenological description of a quasiparticle's 
motion near the NS interface. In Sec. Ill, we discuss the 
conductance peaks in NS junctions. A relation between 
the broken time-reversal symmetry states and the peak 
position in the conductance is discussed in Sec. IV. we 
apply the phenomenological theory to NS junctions un- 
der magnetic fields in Sec. V. In Sec. VI, we summarize 
this paper. 
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II. QUASIPARTICLE'S MOTION NEAR NS 
INTERFACES 



Let us consider two-dimensional NS junctions as shown 
in Fig. ^ where a normal-metal [x < 0) and a super- 
conductor {x > 0) are separated by a potential barrier 
Vif) — VoS{x). We assume the periodic boundary con- 
dition in the y direction and the width of the junction is 
W. The NS junctions are described by the Bogoliubov-de 
Gennes equation;^ 
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FIG. 1: A normal-metal / superconductor junction is illus- 
trated in (a). The trajectories of a quasiparticle in the elec- 
tron branch and those in the hole branch are denoted by solid 
and broken lines, respectively. In (b), the pair potentials of s, 
d^2_y2 and dxy wave symmetries are schematically illustrated. 
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where (ps is a macroscopic phase of the superconductor, 
Rc = {Xc,Yc) = (r + r')/2 and = r - r' . Here we 
assume spin-singlet superconductors for simplicity. The 
argument in the following can be extended to spin-triplet 
superconductors as shown in Appendix. When an elec- 
tronlike quasiparticle is incident from the normal metal 
as shown in Fig. ^ the wave function in the normal metal 
is given by. 
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where and ky are the wave numbers on the Fermi 
surface and they satisfy k'^ + ky = kp with kp being the 
Fermi wave number. Throughout this paper we assume 
that E ~ Aq <^ fj,p^ where Aq is the amplitude of the pair 



potential and E is the energy of a quasiparticle measured 
from the Fermi energy, hf — h'^kp/{2m). In Eq. Q), 



and 



,,he 



are the normal and the Andreev reflection 



coefficients, respectively. 



When a quasiparticle is incident from the normal metal 
in the electron branch, directions of the outgoing waves 
are indicated by arrows as shown in Fig.^ The trajecto- 
ries of a quasiparticle in the electron branch and those in 
the hole branches are denoted by solid and broken lines, 
respectively. In the normal metal, a velocity component 
perpendicular to the interface changes its sign in the nor- 
mal reflection, whereas all velocity components change 
signs in the Andreev reflection. In the superconductor, 
the wave number in the electron branch is (kx,ky), but 
that in the hole branch becomes {—k x^ky^ . In unconven- 
tional superconductors, the pair potential in the electron 
branch (A+ = A(fcx,fcy)) differs from that in the hole 
branch (A_ = A(— fej,, /cy)). Therefore the wave func- 
tion in the superconductor is described by these two pair 
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potentials, 



'J4>- piVs 
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where {t^'^) is the transmission coefficient to the 
electron (hole) branch in superconductors. The wave 
numbers of a quasiparticle are approximately given by 
« fc^ + (-)i/(2^o) fori; - 0, whereto = ?iwF/(7rAo) 
is the coherence length and vp — hkp/m is the Fermi 
velocity. Thus a quasiparticle penetrates into the super- 
conductor within a range of ^o- In Eqs. 101-0, a phase 
g«0± represents the sign (internal phase) of the pair po- 
tential and appears in the wave function in addition to 
a macroscopic phase of the superconductor. The trans- 
mission and the reflection coefficients are obtained from 
the boundary conditions of these wave functions. Near 
the junction interface, an incident quasiparticle suffers 
two kinds of reflection: (i) the normal reflection by the 
barrier potential at the NS interface and (ii) the Andreev 
reflection by the pair potential in the superconductor. In 
this paper, we consider separately contributions of the 
two reflection processes to the reflection coefficients. 

We first consider NS junctions with no barrier potential 
at the interface, 



0, 



(10) 



where zq represents the strength of the potential barrier. 
The Andreev reflection coefficients become 



E-n± 



(11) 
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where 



is the Andreev reflection coefficients from the 



electron branch to the hole branch in the absence of the 
potential barrier. We also give the Andreev reflection 
coefficient from the hole branch to the electron branch 
(rg''). In the case of — Aj_ < 0, v± can be described 
as 



iy± =- 



+ «7 



E 



: cos 9± + i sin 9± = e*^± . 



(14) 
(15) 



Thus the Andreev reflection coefficients include only the 
phase information in the limit of zq = 0. 

We next consider the reflection by the potential bar- 
rier in a phenomenological way. In the presence of the 
potential barrier, the Andreev reflection processes are 
shown in Fig. [2| In the electron branch, the normal 
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FIG. 2: The Andreev reflection processes are decomposed 
into a series of the reflections by the pair potential and the 
barrier potential. 

transmission and the normal reflection coefficients of the 
barrier are calculated to be In = kx/{kx + izo) and 
rjv = —izo/{kx + izo), respectively with k^ = k^/kp- 
Those in the hole branch are t*^ and r'^j. The Andreev 



reflection coefficient in the 1st order process is given by, 



(16) 



At flrst an electronlike quasiparticle starting from 
transmits into the superconductor through rp, (ijv)- In 
Fig- 13 the vectors in real space are surrounded by squares 
to avoid confusion. While traveling the superconductor 
within the range of ^q, the quasiparticle is reflected into 
the hole branch by the pair potential at ri {r^^)- Then 
the quasiparticle goes back to the normal metal in the 
hole branch through r'^ (i^). The 2nd order Andreev 
reflection process in Fig. [3 (b) can be estimated in the 
same way, 
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After the first Andreev reflection into the hole branch, the may suffer additional phase shift while moving around 
quasiparticle suffers the normal reflection (r^). Next the the NS interface. Actually, an electron acquires a phase 
holelike quasiparticle experiences the 2nd Andreev reflec- ^ik-in-ro) -^j^jig traveling from rp to ri as shown in Fig. [3 
tion to the electron branch at r2 Then the elec- (a). In addition to this, a phase factor e*'" '-''''"''^^ is mul- 

tronlike quasiparticle suffers the normal reflection (vn) 
followed by the 3rd Andreev reflection into the hole 
branch {tq'^)- Finally the holelike quasiparticle goes back 
to the normal metal through Tq (t^). We only show the 
expression of the Andreev reflection coefficient in the 3rd 
order process. 



= (3) 



+* .42 he 



t]y. 



(20) 



The corresponding trajectory is shown in Fig. [21(c). The 
total Andreev reflection coefficient is obtained by the 
summation of these reflection processes up to the infl- 
nite order. 



tiplicd while traveling from ri to Tq in the hole branch. 
These two phase factors exactly cancel each other out 
when the retro property holds because Tq = Tq. Thus 
rf^ indicates the same position for all n. In particu- 
lar for E = 0, a, relation r.^ — r for all n holds, which 
means the retro property of a quasiparticle in the normal 
metal. In the limit of E 0, we find in Eq. H15|l that 
I'i — > 1 irrespective of symmetries of the pair potential. 
The Andreev reflection probability becomes 
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In the similar way, the normal reflection coefficient results 
in 



Firstly, we consider superconductors where the pair po- 
tentials in the two branches (A+ and A_) have the same 
sign, (i.e., e*^'^"^'^+^ = 1). For examples, the pair poten- 
tials below satisfy the condition irrespective of the wave 
numbers of a quasiparticle, 
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Ad ,(fc) -Ao(fc2 



(4 



(s wave), (26) 
wave), (27) 



Although the reflection coefficients in Eqs. H21() and H22() 
are obtained based on the phenomenological description 
of a quasiparticle's motion, they are mathematically iden- 
tical to the exact expressions calculated from the bound- 
ary conditions of the wave functions in the presence of 
the potential barrier^ 



III. CONDUCTANCE 

The differential conductance is calculated from the nor- 
mal and the Andreev reflection coefficientsi2i*2^ 
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E=eVti. 



where Vtias is the bias voltage applied to NS junctions. 
We focus on the limit of i? — > for a while, where the 
Andreev reflection probability dominates the zero-bias 
conductance because the conductance can be described 
by 



= ^ E I 
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E=eVti 



A quasiparticle after the Andreev reflection traces back 
the original trajectory of a quasiparticle before the An- 
dreev reflection. This is called the retro property of a 
quasiparticle. When we estimate the reflection coeffi- 
cients in Eq. H21|) and (|22|l . we only consider the phase fac- 
tor of the Andreev reflection. A quasiparticle, however. 



where — kx/kp and ky — ky/kp are the normalized 
wave number on the Fermi surface in the x and y di- 
rections, respectively. The schematic flgures of the pair 
potentials are shown in Fig. ^ (b) . Equation H2t)|) repre- 
sents the pair potential of s wave superconductors. The 
pair potential in Eq. (|27|l is realized in a junction where 
the a axis of a high- Tc superconductor is perpendicular 
to the interface normal. When e*^'^"^'^+'' = 1 is satis- 
fled, Eq. H25() becomes the summation of the alternating 
series. The Andreev reflection probability results in 
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In low transparent junctions, (i.e., 3> 1), the Andreev 
reflection probability becomes a small value |iAr|*/2 cx 
l/^Q. Therefore the zero-bias conductance in Eq. H24|l 
is proportional to l/zg. Secondly we consider that the 
signs of the two pair potential are opposite to each other. 
The pair potential 



(fc) 



2AQk,.ky, 



{dxy wave). 



(29) 



satisfies 6''-'^""'^+-' = — 1 for all wave numbers and is real- 
ized in a junction where the a axis of a high- Tc supercon- 
ductor is oriented by 45 degrees from the interface nor- 
mal. All the expansion series in Eq. (|25|l have the same 
sign and the Andreev reflection probability becomes 
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1. 



(30) 



Thus the zero-bias conductance in Eq. (|24|) takes its max- 
imum value. The sign of the pair potentials character- 
izes the interference effect of a quasiparticle near the NS 
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interface. For 6*'-'^""'^+^ = 1, the alternating series in 
Eq. H25I) reflect the destructive interference among the 
partial waves of a quasiparticle in the expansion series. 
Hence the conductance becomes small at the zero-bias. 
On the other hand for e*^'^-~'*+-' = —1, the expansion 
series with the same sign imply that the partial waves in 
the expansion series interfere constructively, which leads 
to the large zero-bias conductance. The constructive in- 
terference at the interface causes a resonant state which is 
now referred to as the ZES. The Andreev reflection prob- 
ability is unity independent of the normal transmission 
probability of junctions as shown in Eq. H3U|I . This can 
be interpreted as a result of the resonant transmission of 
a quasiparticle through the ZES. A microscopic calcula- 
tion shows that the ZES has a large local density of states 
around x = (^q at the zero-energy.^^ Similar arguments 
have been done in normal- metal/insulator/normal- 
metal/insulator/superconductor junctions^'* and at the 
surface of high-T^ superconductors.^^ 

In Eq. (|30|) . we can explain a large conductance at 
the zero-bias. In what follows, we will show that the 
conductance has a peak structure around the zero-bias. 
When E but still E ^ Aq, the degree of resonance 
is suppressed because i/± is no longer unity as shown in 
Eq. IjlSfl . In the superconductor, the argument of the 
phase cancellation in the round-trip between Tq and ri 
in Fig. 121(a) is still valid as far as E^ — |A±p < being 
satisfied. In the electron branch on the way to ri, the x 
component of the wavenumber is given by 



(31) 



kx "^fJ-F 

The real part determines the direction of the quasiparti- 
cle's motion. The inverse of the imaginary part charac- 
terizes the dumping of the wave function and is roughly 
estimated to be ^o- It is also shown that kx is the real 
part of the wave number in the hole branch on the way 
back to tq. The Andreev reflection probability for finite 
E is given by 

\tN\' 



„he\2 



(32) 



To make clear a relation between the peak positions of 
the conductance and the relative sign of the pair poten- 
tials, we consider the pair potential 

^rf.iE„(fe) = ^0 Sgn(fc::cfcy), (33) 

instead of Eq. 1)29(1 . Here the anisotropy of pair potential 
is taken into account only through the phase e*"^* and 
the k dependence of the pair potential is neglected. The 
pair potential in Eq. (|33|l is illustrated in Fig.n](b). We 
will check the validity of Eq. 1(33(1 later. The Andreev 
reflection probability for Ad^j „ becomes 
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(35) 



where we use a relation O = Oj^=0-\n Eq. 1(15(1 . The 
Andreev reflection probability has a peak structure at 
E = Q and the width of the peak is characterized by Eq 
which is Ao/zq in the limit of Zg ^ 1. On the other hand 
in s wave junctions (i.e., e*^'^-~'^+^ = 1), we find 
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The Andreev reflection probability has a peak at £' = Aq 
reflecting a peak of the bulk density of states in s wave 
superconductors. In Fig. |31 we plot the conductance, 
where zq = 3 and Nc = Wkp/n is the number of prop- 
agating channels on the Fermi surface. The results for 
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FIG. 3: The conductance is plotted as a function of E, where 
zo — 3. The anisotropy of the pair potential in d^y symmetry 
is taken into account only through the phase factor e'^'^~~'^+' 
and k dependence of the pair potential is neglected in the 
dash-dotted line. The d^j/ symmetry is fully taken into ac- 
count in the solid line. The conductance for d^i_y2 symmetry 
is amplified by 5 times in the dotted line. 

s wave junction are indicated by the broken line. The 
conductance for dj.2_y2 symmetry in the dotted line is 
amplifled by 5 times. These conductance has a peak at 
E = Aq reflecting the bulk density of states. The re- 
sults for Ad^j^jj and dxy are shown with the dash-dotted 
line and the solid line, respectively. There is no signifl- 
cant difference between the conductance for d^y symme- 
try and that for dgign because the relative sign of the two 
pair potential (^e^('f'-^ '!'+'> — —1) dominates the subgap 
conductance structure. Throughout this paper, we de- 
scribe the pair potential by using the step function at 
the NS interface and neglect its spatial dependence in 
superconductors. In real NS junctions, the pair poten- 
tial is suppressed at the interface in the presence of the 
2jT;g,60j65 r^Yie conductance shape around the zero-bias, 
however, almost remains unchanged even if the spatial 
dependence of the pair potential is taken into account 
This is also because relative sign of the two pair poten- 
tials determines the conductance around the zero-bias. 
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The spatial dependence of the pair potential may affect 
the width of the ZBCP through Eq in Eq. 

We note that there is no remarkable differences be- 
tween the mathematical origin of the peaks at E ^ for 
gi(0_-</,+) ^ _i j^jjjj ^.jjj^i- at £; = Ao for e*('^--'^+) = 1. 

Actually it is easy to confirm at E — Aq that the Andreev 
reflection probability in s wave junctions becomes 
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n=0 



Ji<P~-4>+) 



(37) 



All the expansion series have the same sign for 

gi(0--0+) = 1^ 

In above arguments, we have assumed that the junc- 
tions have non zero transmission probabilities. In the 
end of this section, we briefly mention that the ZES be- 
comes a real bound state in the limit of 2:0 — > oo. A 
quasiparticle motion is spatially limited at the surface 
of the semi-finite superconductor because of the perfect 
normal reflection by the surface and the Andreev reflec- 
tion by the pair potential. The ZES becomes a bound 
state because there is no quasiparticle excitations which 
extend into the bulk superconductors at £^ = 0. In the 
density of states, such ZES is found as the (5-function 
peak. For finite transmission probability of junctions, 
the finite propagation into normal metals gives a finite 
life time of the ZES which is given hy h/Eg. On the other 
hand for e*^'^-^'*+) = 1, the resonant state at E — Aq 
does not become a bound state because there are excita- 
tions extend into the bulk superconductors at E = Aq. 
In superconductor /insulator /superconductor (SIS) junc- 
tions, the ZES is also a bound state irrespective of the 
transmission probability of junctions. The description of 
the Andreev bound states in SIS junctions was given, for 
example, in Ref. IttI 



with a^ + (3^ = 1. We find 



|A±| 

-0+) 



a^Al 



(39) 
(40) 



In Fig ^ we show the conductance in the s + idxy sym- 
metry for several a. For a = 0, the results are identical 
to the conductance of dxy wave junctions in Fig. |31 The 
ZBCP splits into two peaks for a 7^ 0. The splitting 
width increases almost linearly with increasing a. In the 
limit of a = 1, the results coincide with the conductance 
of s wave junctions in Fig. |21 The peak position can 
be explained by the expression of the Andreev refiection 
probability 
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|iw|4 + 4|rw|2cos2(0 + ,^_)' 



(41) 



/I A 12 _ jp2 p 

cos{9 + </._) aAo + ^/?A,,„, (42) 



« + 7^/3 sgn(fc,fc,). (43) 

In the last equation, we replace A^^^^ by A^^^^^. The 
conductance peak (the resonance energy) is expected at 
an energy which satisfies cos{9 + <j}-) — Q as shown in 
Eq. (|41|l . The resonance energies for a — Q and a — 1 
are E = Aq and E — respectively. These resonance 
energies are independent of the wave numbers. Conse- 
quently the peak heights for a = and a = 1 become 
unity. The peak heights for finite a, however, are always 
less than unity as shown in Fig.0]because the resonance 
energy depends on wave numbers as shown in Eq. H42|l . 



IV. PAIRING WITHOUT TIME-REVERSAL 
SYMMETRY 

In recent experiments, a possibility of the broken time 
reversal symmetry state (BTRSS) at the surface of high- 
Tc superconductors has been discussed i^^i^SiSi^SiS These 
experiments found the split of the ZBCP in the zero mag- 
netic field. It is pointed out that such surface states may 
have s-\-idx^^ or dxy + id^^-y"^ pairing symmetry. The- 
oretical studies showed the split of the surface density of 
statea^^iS24SiSliSiS& when s+idxy wave pairing symmetry 
is assumed at the surface of the dxy wave superconduc- 
tor. Within the present phenomenological theory, it is 
also possible to discuss the split of the conductance peak 
by the BTRSS in terms of the shift of the resonance en- 
ergy. We assume the pair potential as 

As+id^^{k) = aAo + if] Ad^^{k), {s + idxy wave) , (38) 




-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 

E/A, 



FIG. 4: The conductance is plotted as a function of E for 
s + idxy symmetry, where zo — 3. 

The positions of the conductance peaks are roughly 
given by _B = ±aAo, which can be understood by the 
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resonance condition of cos{6 + (p-) = in Eq. H43|) . Since 
the peak position is determined by a, relative ampli- 
tudes of s and d^y components can be estimated from 
the peak splitting width observed in experiment. In the 
phenomenological theory, effects of the BTRSS on the 
conductance can be understood in terms of the shift of 
the resonance energy. 

In theoretical studies, it is shown that the s + id^y 
wave BTRSS splits zero-energy peak of the local den- 
sity of states^MSiSiSM^ and the ZBCPS Experimen- 
tal results are, however, still controversial. Some ex- 
periments reported the split of the ZBCP at the zero 
magnetic field ,SiS&SLSSiS other did not observe the 
splittingi''i^°i^^i^'^i-^^i^''i^^ Thus opinions are still divided 
among scientists on the BTRSS in high- Tc superconduc- 
tors. If the BTRSS does not exist, we have to find 
another reasons for the peak splitting observed in ex- 
periments. In recent papers, we have showed that the 
interfacial randomness causes the split of the ZBCP in 
the zero magnetic field in both numericalljiS using the 
recursive Green function methodSiSS and analyticallj^ 
using the single-site approximation^ Our conclusion, 
however, contradicts to those of a number of theories 
^2i22i2MSi2SiSi based on the quasiclassical Green function 
method.2Si22i2ai2ii2^ The drastic suppression of the ZBCP 
by the interfacial randomness is the common conclusion 
of all the theories. The theories of the quasiclassical 
Green function method, however, concluded that the ran- 
dom potentials do not split the ZBCP. 



V. EFFECTS OF MAGNETIC FIELD 

The TRS is also broken by applying external magnetic 
fields onto NS junctions. The resonance at i? = is 
suppressed because a quasiparticle acquires a Aharonov- 
Bohm like phase from magnetic fieldsiSiSiS^ Actually it 
is pointed out that the ZBCP in NS junctions splits into 
two peaks under the magnetic field The re- 



N 



N 



e 1-1= (x y) 




r2=(x2,y2) 




(a) 



(b) 



FIG. 5: The motion of a quasiparticle near the interface is 
illustrated. 

fiection process in Fig.|S{a) corresponds to As in Eq. (|19|l . 



We consider uniform magnetic fields perpendicular to 
the xy plane (i.e., Bz) and assume the Landau gauge 
Aext = Bxy. Effects of magnetic fields are taken into 
account through the phase of the wave function by using 
the gauge transformation. While traveling from rp to ri, 
an electronlike quasiparticle acquires a phase e*"^™ with 



e 

he 



ri 

dl-AextH) 



e B 
hc'2 



{xi+xo){yi~yo). (44) 



Since the magnetic field is sufficiently weak, the integra- 
tion path can be replaced by a straight line between rp 
and ri which is denoted by Ci in Fig. I3a), This ap- 
proximation is justified when the radius of the cyclotron 
motion of a quasiparticle, 2iJ,p/{kpheB /mc), is much 
larger than ^o- The condition is equivalent to the relation 
ttAq ^ heB/mc. In high- Tc materials, Ag ~ 30 — 40 
meV, whereas heB/mc is 10"^ meV for B = 1 Tesla, 
where we use the bare mass of an electron. The phase 
shift on the way from ri to Tq (C2) in the hole branch is 
equal to e"^™. This is because the direction of a quasi- 
particle's motion and the sign of the charge on Ci are 
opposite to those on C2 at the same time. In the same 
way, we can show that the phase shifts on C3 and C4 in 
Fig. |31 (a) are also e*"^™ . Under the gauge transformation, 
the pair potential should be changed to 



A(r, r') exp 



At ri, a phase factor 




dl ■ Aext{l) 



exp 



~i2e 

he 



dl ■ A^xt{l) 



(45) 



(46) 



is mulitiplied to the Andreev reflection coefficients, where 
r and r' in Eq. H45I) are set to be ri. A phase factor 



exp 



i2e 
he 



r-2 

dl ■ Aext{l) 



(47) 



is also multiplied to the Andreev reflection coefficients 
at r2- The total phase shift by the magnetic field along 



C4 in Fig. I^Ia) (e*^"^^ ) is then given by 



e 

he 



r2 

dl ■ Aext{l), 

ri 



eB 

he 



[yi ~ yo)(a;i - a;o) 



B ky 
Bo kx ' 



So 



(48) 
(49) 
(50) 



where (po = 2Tihe/e. On the way to Eq. I|49|l . we use a 
relation {xi - xo)/{yi - yo) = kx/ky and xi - xq ^ ^q. 
We note that 2(j)B is the gauge invariant magnetic flux 
passing through the gray area in Fig. [Sfb) , where — 
(2xi + Xo,yo). Thus 2(f)B remains unchanged in another 
gauges such as A^xt — ~Byx and penetrating magnetic 
flelds Aext = BXoe~^^^°y with Aq ^ C07 where Aq is the 
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FIG. 6: The conductance under magnetic fields for d^y sym- 
metry, where zo = 3 and 10 in (a) and (b), respectively. In 
high-Tc material, Bo is estimated to be 160 Tesla. 



penetration depth. In high- Tc materials, ~ 2nm and 
Ao 200nm. 

Effects of magnetic field can be taken into account in 
the present theory by 



(51) 



where Ag is defined in Eq. (|19|l . We show the conduc- 
tance in dxy wave junctions calculated from Eqs. pi() - 
(P^ and in Fig. El where zq = i and 10 in (a) 
and (b), respectively. In high- Tc superconductors, Bq is 
about 160 Tesla. The ZBCP decreases with increasing B 
in both Figs.El(a) and (b). The degree of suppression due 
to magnetic fields depends on the transmission probabil- 
ity of the junction. More drastic suppression can be seen 
in lower transparent junctions. In Fig. El (b), the ZBCP 
almost disappears for B = 0.05i3o- The ZBCP, however, 
remains one peak and does not split into two peaks even 
in the strong magnetic fields. The results in Fig. El are 
qualitatively well described by the analytical expression 



of the Andreev reflection probability for ^ Aq , 



(52) 
(53) 



|tAr|4 + 4|rjvPsin2(0 + 0B)' 

'\tM\W + '^\rN?{E+\^\^Bf 



We linearize the magnetic fields in sin(6'-f 0b) in Eq. H53|l. 
Eauation (|53(l implies that the resonance energy may be 
shifted from i? = by magnetic fields. In contrast to the 
sphtting of the ZBCP by the BTRSS in Sec. IV, we do 
not find the peak splitting under magnetic fields in Fig. El 
In the BTRSS, the shift of the resonance energy is caused 
by the s wave component which has the resonance energy 
at E = Aq. On the other hand, any resonant states are 
not associated with magnetic fields. Thus the magnetic 
fields only suppress the resonance of the ZES as shown 
in Fig. El 

In a previous paperj^'* however, the split of the ZBCP 
in magnetic fields was reported within the quasiclassical 
approximation (QCA). The results in Eq. (|53(l are sim- 
ilar to that in the argument of the Dopplar shift in the 
QCA. The supercurrents flows along the interface shift 
the energy of a quasiparticle as 



c 



Ps = = -e-^/^"y, 



(54) 
(55) 



where is the condensate momentum at the interface. 
In Eq. H55(l . d wave character of the supercurrent is 
not considered. The corresponding approximation in the 
present theory is replacing E + \ A\(j)B hy E + Aqc^b in 
Eq. fK!!| and we find 



„he\2 



\tNm 



I QCA 



|i^|4|A|2+4|rjvP(i? + Ao</)B)2 



(56) 



In Fig. [71 (a), we show the conductance calculated from 
Eq. for zq = 10. In contrast to Fig. El (b), we find 
split of the ZBCP when magnetic fields are larger than 
the threshold magnetic field, Be- The threshold depends 
on zq as shown in Fig. [7\ (b), where Be is plotted as 
a function of 1/zq which is proportional to the normal 
transmission probability of junctions. The threshold in- 
creases with increasing the transmission probability of 
junctions. This has been pointed out in the conductance 
calculated on the lattice model by using the QCAi^ In 
the lattice model, it was also shown that Be decreases 
with the increase of the doping rate. The Fermi energy 
is a decreasing function of the doping rate. Therefore 
the transmission probability of junctions decreases with 
increasing the doping rate. 

Although Eqs. (|53|l and H56|l are similar to each other, 
the response of the ZBCP to magnetic fields are qualita- 
tively different. To make clear if a magnetic field splits 
the ZBCP or not, we need some numerical simulations, 
where effects of magnetic field are taken into account ac- 
curately. In experiments, some papers show the split of 
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VI. CONCLUSION 




-0.05 0.00 0.05 



E/A„ 



FIG. 7: The conductance in the quasiclassical approximation 
is plotted for several magnetic fields in (a), where zq — 10. 
In (b), threshold magnetic fields are shown as a function of 
1/zl 



We have presented a phenomenological theory of the 
Andreev reflection to make clear reasons for the ap- 
pearance of the zero-bias conductance peak (ZBCP) 
in normal-metal / unconventional superconductor junc- 
tions. The phenomenological theory reveals that the 
zero-energy state (ZES) is a consequence of the construc- 
tive interference effect of a quasiparticle. The expression 
of the Andreev reflection probability enables us to under- 
stand an importance of the unconventional pairing sym- 
metry for the formation of the ZES. The phenomenolog- 
ical theory is applied to superconductors with a broken 
time-reversal symmetry state (BTRSS) and junctions un- 
der magnetic fields. The split of the ZBCP in s + id^y 
wave superconductors is understood in terms of the shift 
of the resonance energy by the s wave component. The 
Aharonov-Bohm like phase received from magnetic fields 
suppresses the degree of resonance of the ZES, which ex- 
plains the suppression of the ZBCP in magnetic fields. 



APPENDIX A: ANDREEV REFLECTION BY 
SPIN-TRIPLET SUPERCONDUCTORS 



the ZBCP in magnetic fieldsi^^*^ On the other hand, sev- 
eral papers report no splitting of the ZBCP.iSi2Si2L2a A 
microscopic scattering theory indicates that the sensitiv- 
ity of the ZBCP to magnetic fields depends on the degree 
of potential disorder near the NS interfaceiS 



In the text, we consider two-dimensional spin-singlet 
superconductors and 5-function type potential barrier for 
simplicity. Here we generalize the phenomenological the- 
ory to spin-triplet superconductors in three-dimension. 
The pair potential in superconductors is given by 



Finally we briefly discuss an important difference of 
the conductance in the present theory and that in the 
QCA. The phenomenological theory reaches at Eq. H56() 
which is almost the same as the conductance expression 
in the QCA.^"' The two theories, however, still lack a 
quantitative agreement of the threshold magnetic field. 
The normalization for the penetrating magnetic fields in 
the QCA (B^'^^ = 0o/(2<oAo)) is about 1.6 Tesla with 
Ao - 100^0-^^ In Eq. ^,p^m the QCA is originally 
given by the vector potential which is not an observ- 
able value. Thus the QCA does not satisfy the gauge 
invariance. In the present theory, on the other hand, 
we consider uniform magnetic field and the normaliza- 
tion of magnetic fields (Bq) is about 160 Tesla. This 
value remains unchanged even if we consider penetrating 
magnetic field as Be~^^^°z with Aq ^ £,o- For exam- 
ple in Fig. (|7J) (b), we find that Be is about O.IBq at 
Zq = 3. Therefore Be is estimated to be 16 Tesla in 
the present theory. The same results are interpreted as 
Sc=0.16 Tesla ii we use B^^"^ in the QCA. The thresh- 
old magnetic field in the QCA is estimated to be much 
smaller than that in the present theory. This disagree- 
ment may be important because the maximum value of 
magnetic fields in experiments is about 10 TeslaiiSi2S42L2& 



A(fe) = 



id{k) ■ a-&2 ■ triplet. 



id{k)a2 



singlet. 



(Al) 



where aj for j = 1, 2 and 3 are Pauli matrices repre- 
senting the spin degree of freedom. We assume that the 
current is in the x direction and consider a potential bar- 
rier 



V{r) = Vo [e(x) - e(.T - L)] 



(A2) 



where L is the thickness of the insulating layer. The An- 
dreev reflection coefficients in the absence of the insulator 
are calculated analytically 



(-)^(-)' 



A(±) —id± ■ o-(T2, 



L± 



A2 



Pi 



(A3) 
(A4) 
(A5) 

(A6) 
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Pi,± =|q±|o-o - ("1)'9± ■ d-, 
=id± X d^, 
d± =d{±k.^,ky,kz), 



(A7) 
(A8) 

(A9) 
(AlO) 
(All) 



where ips is a macroscopic phase of superconductor, l{= 1 
or 2) indicates the two spin branches of Cooper pairs and 
(To is the 2x2 unit matrix. The normal transmission 
and the normal reflection coefficients of the insulator are 
calculated as 



In 



ik^ L 



14) 



(A12) 
(A13) 
(A14) 



zq = sinh(p^L), 

^1 ^{pI - K:) sinh(p^L) + 2ikxpx cosh(p^L), (A15) 



where Px = \/ {Vq/ ^j.f) — (kx/kp)^ is the wave number 
at the insulator and px = Px/kp- 

The argument in Sec. II leads to the exact expression 
of the Andreev and the normal reflection coefficients^^ 
which are given by 



r = - zqZi 



r = — e 



\z^\^o^-zlW 



Aklpl^\^.R(+) \z^\^o^~zlW 



W^=^(-)A(-)A[+)i?(+). 



(A16) 

-1 

7 

(A17) 
(A18) 



The results of unitary states including the spin-singlet 
states can be obtained when we use following relations 



R 



(±) 



E 



\D±\ 



in, I _ / M±l ■ singlet 
|d±| : triplet, 



(A19) 
(A20) 



in Eqs. ljA16|) - (jA18|) . The differential conductance is 
given by 



G 



NS 



h 



/le/-;/le^t^ 



(A21) 



A relation d_ = — d+ represents the condition for the 
perfect formation of the ZES. Actually when d+ — d — 
I'd- with v = ±1, the Andreev reflection probability 
becomes 



Tr [f^'^if''^)^] = 



1=1 



ikiPiAiKi 



AklPlAf + ziiAf-jyKf 



where Ki ~ Ki,+ = Ki.- and A; = A/.+ = A/_, 
In the limit of ^ and zq ^ 1, we find 



(A22) 



Tr [f'^^'if'^y] = 



(A23) 



where spin degree of freedom give rise to a factor 2. Thus 
the zero-bias conductance is independent of the transmis- 
sion probability of junctions when d = d^ is satisfied. 

In spin-singlet superconductors, we show that the in- 
ternal phase of a Cooper pair is responsible for the ZES. 
In spin-triplet superconductors, the internal spin degree 
of freedom of a Cooper pair has another possibilities for 
the formation of some resonant states in subgap energies. 
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